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We study the equilibrium dynamics of the relative phase in a superconducting Josephson link 
taking into account the quantum fluctuations of the electromagnetic vacuum. The photons act as a 
superohmic heat bath on the relative Cooper pair number and thus, indirectly, on the macroscopic 
phase difference (j). This leads to an enhancement of the mean square {(j)'^) that adds to the spread 
due to the Coulomb interaction carried by the longitudinal electromagnetic field. We also include the 
interaction with the electronic degrees of freedom due to quasiparticle tunneling, which couple to the 
phase and only indirectly to the particle number. The simultaneous inclusion of both the radiation 
field fluctuations and quasiparticle tunneling leads to a novel type of particle-bath Hamiltonian in 
which the quantum particle couples through its position and momentum to two independent bosonic 
heat baths. We study the interplay between the two mechanisms in the present context and find 
interference contributions to the quantum fluctuations of the phase. We explore the observability 
of the QED effects discussed here. 

I. INTRODUCTION 

A proper understanding of decoherence as resulting from the dissipative effect on a quantum system of a complex 
environment is of paramount importance in diverse areas of modern physics and technology. The influence of many 
degrees of freedom on a given quantum variable has for long been considered a problem of fundamental interesti, 
partly because of its connection to the polaron problem^. However, it was only in the eighties when, because of its 
relevance to macroscopic quantum tunneling and coherence^, the problem was addressed by a large community of 
physicists. The loss of quantum wave coherence -also referred to as decoherence or dephasing- appears in a wide range 
of physical contexts. It has been generally recognized that decoherence acts much faster than energy dissipation, since 
it only requires the excitation or destruction of a single quantum of the thermal oscillator bath, while many quanta 
are necessary to change the particle energy appreciablji^iSii. In the context of mesoscopic electron systems, dephasing 
is responsible for the destruction of those quantum interference effects which characterize transport through ballistic 
nanostructures or disordered metals at low temperatures. Particle-bath models have been used to address such a 
problemSi^. In disordered metals at sufficiently low temperatures, the dominant mechanism responsible for the loss of 
phase coherence of electrons near the Fermi surface is the Coulomb interaction with identical electrons in the metalSi. 
The question of whether or not the electron dephasing time diverges at zero temperature has been the subject of a 
recent controversyifi. 

We wish to emphasize here that, under the generic term decoherence or dephasing, one may be contemplating 
physical problems of widely different conceptual and technical nature which may yield different answers to apparently 
similar questions. For instance, the general question 'does zero-point decoherence exist?' does not admit a universally 
valid response. The answer depends on the system under attention, the nature of its surrounding environment, the 
properties of the coupling, and the physical representation within which we frame the question. The latter remark 
refers to the importance of the basis set in whose representation the reduced density matrix of the quantum system 
tends to become diagonal as a result of the interaction with the environment. The choice of representation is related 
in turn to the observable to be measured, since the expectation value must be sensitive to the off-diagonal terms 
of the density matrix in the desired representation. As examples of distinct physical problems having to do with 
zero-point fluctuations, we may mention the problem of decoherence for a quantum particle coupled to a dissipative 
bath, which is qualitatively different from that posed by Fermi surface electrons interacting with other electrons in 
a disordered metal. An intermediate class of problems is that of conduction electrons interacting with an external 
environment, which has been shown to modify some Fermi liquid properties such as persistent currents in mesoscopic 
ring&ii. The message is that possible answers to questions having to do with the loss of quantum wave coherence 
must be investigated case by case and cannot be translated from one physical context to another without a careful 
scrutiny. 

The general problem of decoherence is receiving renewed attention due to its central role in the design of quantum 
information processing systems. In this context, it seems that all potential sources of decoherence must be explored. 
Due to its ubiquity, the quantum electrodynamic (QED) field provides the most basic decoherence mechanism that a 
charged particle may experience. The purpose of this paper is to investigate the role of the quantum electromagnetic 
(EM) field as a possible source of decoherence. Specifically, we aim at understanding the effect of the transverse 
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component of the EM field. It must be distinguished from the longitudinal component, which in the radiation gauge is 
originated by the electrostatic Coulomb interaction among charged particlesiSi. The combined effect of the transverse 
and the longitudinal field in a many-body context has been recognized to be responsible for logarithmic departures 
from the Fermi liquid picture in two-dimensional electron systems^^ and for the reduction of persistent currents in 
mesoscopic rings^*. In the context of particle-bath problems, it is fair to say however that, despite its fundamental 
character, the QED field has been little explored as a dissipative environment that may disturb electrons and, in 
general, charged particles. From the quantum dissipation viewpoint, Quantum Electrodynamics poses at least two 
new problems. First, in vacuum it is a superohmic environment^^'^^ and, as such, it has been much less studied than 
its ohmic counterparts. Second, it has no intrinsic upper cutoff. Thus, when needed, appropriate cutoffs must be 
introduced on the basis of sound physical arguments. 

Some studies can be found in the literature dealing with the coherence properties of free charged particles in the 
presence of the background radiation fieldi^4i2*i2iiS. To understand in greater depth the effect of the transverse photon 
field as a source of dissipation, we choose as our case study a system which is otherwise quite well understood: the 
phase-number variable of a Josephson link connecting two superconductors. Typically, such a weak link is achieved 
with a tunneling barrier or a point contact. The problem which we wish to address here may be viewed as the 
macroscopic quantum-mechanical version of the Lamb shift problem, since we aim at calculating the effect of the 
fluctuations of the transverse photon field on an otherwise conservative quantum-mechanical problem. In particular 
we will focus on the calculation of the uncertainty At/) = {(fy^) — (0)^, where cj) is the relative phase between the two 
superconductors. We will focus on the equilibrium dynamics, for which (^) = 0, and on the harmonic limit {(j)'^) <t:l. 
The effect of the transverse EM field on {(j)"^) can be potentially detected, for instance, as a Debye- Waller reduction of 
the Josephson critical current or as an intrinsic source of uncertainty in the fine measurement of tiny magnetic fields. 

This paper builds on a preliminary study presented in Ref^, where it was shown that, in the absence of a Josephson 
coupling, recurrent dynamics alternating between quantum collapses and revivals is robust against fluctuations of the 
QED field^ii. Such collapses and revivals were later predicted^^ and observed^^ in atomic Bose-Einstein condensates. 
However, a number of theoretical questions, including the role of the upper frequency cutoff, were left unexplored. 
In the present work, we also investigate the effect of quasiparticles, which cause fluctuations of the longitudinal field 
and which we approximate by effective oscillators, adapting the work of Ref.^'^ to a Hamiltonian description. When 
computing the combined effect of both the QED and quasiparticle fields, we find that, in the denominators of the 
spectral weight contributing to (c/)'^), there are terms which are caused by the interference between the two baths. 

In Ref.'^''' current and voltage fluctuations in tunnel junctions were studied. There special attention was paid to the 
biased case and to the effect of fluctuations on the I — V characteristics. Here we have focussed on the calculation 
of the equilibrium values of {(jP) and {N'^) and have used a consistent particle-bath approach. Moreover, we have 
attempted to clearly separate the transverse and longitudinal electromagnetic fields and have modelled quasiparticles 
with effective oscillator baths. Our work presents a compact and somewhat involved treatment of the effect of the 
photon and quasiparticle fields on the phase of a superconducting weak link. Fortunately, some aspects of this complex 
problem may be understood in simpler terms for some particular cases. This observation motivates section II, which 
is devoted to provide a largely self-contained preview of some of the results to be derived more rigorously later. That 
discussion is inspired to some extent in simplified derivations of the atomic Lamb shift (see e.g. RefsSS,). Starting from 
the classical Langevin equation satisfied by 0, the EM field is introduced by invoking gauge invariance arguments. 
In the simplest cases, we are able to derive fourth-order differential equations which, in the absence of quasiparticles, 
reduce to Abraham-Lorentz equations for both N and (f>. We argue conclusively that the coupling to the EM field is 
essentially that of a fluctuating dipole. 

Section HI deals in detail with the coupling of the macroscopic superconducting phase to the QED field, which is 
obtained invoking gauge invariance and introducing appropriate canonical transformations. The number and phase 
autocorrelation functions are calculated. This permits to compute the QED correction to the longitudinal field resultS& 
for (0^) in terms of the junction aspect ratio and the fine structure constant. In section IV, we study the coupling to 
the quasiparticle bath, i.e. to the many-body fluctuations of the longitudinal field. We adapt the work of Ref.^'^ to 
a Hamiltonian language where the dynamics is investigated through the Heisenberg equations of motion. Section V 
addresses the most complex issue, namely, the combined effect of the photon and quasiparticle fields. There we face 
a fundamental problem in quantum dissipation that has so far received little attention: the behavior of a quantum 
particle interacting with two different baths through its position and momentum^^. Finally, section VI is devoted to 
a concluding discussion. 

II. GENERAL REMARKS ON THE EQUATIONS OF MOTION 

The standard way to describe the dynamics of the macroscopic phase in Josephson links is the resistively and 
capacitively shunted junction (RCSJ) model2S*22iS2iffi, which contemplates an ideal Josephson junction shunted by a 
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resistor and a capacitor. The resistor models the dissipative effect of incoherent quasiparticle tunneUng through the 
junction, while the capacitor accounts for the charging energy, which plays the role of a kinetic energy for the phase. 
In the absence of driving currents, the RSCJ model reads 

N{t) = -^sinm^^>p{t) + ^J{t) (2.1) 

0(0 = ^Nit), (2.2) 

where N, the number of transferred Cooper pairs across the junction, is canonically conjugate to the relative phase 0. 
We have introduced the notation Ej = Tilc/ie and Ec = 4e^/C for the Josephson coupling energy and the charging 
energy respectively {Ic and C are the critical current and the capacitance of the junction). /(<) is a stochastic process 
with zero mean. At high temperatures it is related to the resistance by Einstein's relation (/(i)/(0)) — 2kBT6{t)/ R, 
where R is the resistance of the junction in the normal state. Eq. H2.2|) is recognized as the ac Josephson relation 
written in terms of the relative Cooper pair number, which generates a chemical potential difference through the 
interaction energy Ec = d^/dN. 

The effect of the EM vacuum fluctuations may be introduced through the following argumentation: Written in 
language of the Coulomb gauge (as is standard in these physical contextsSSiSS), Eqs. (|2.1I) and (|2.2I) relate the phase 
to gauge invariant quantities. More specifically, they include the effect of the longitudinal electric field, which in 
its simplest form yields a circulation E|| • dr = 2eN/C. Therefore, we are entitled to replace the phase in l|2.2|l 

by its gauge invariant expression cj) = — (j)2 + (2e/hc) dr • A(r), where 1 and 2 are points deep enough in the 
bulk of superconductors 1 (left) and 2 (right). The resulting equation is again interpreted in the particular language 
of the transverse gauge, in which the vector potential is characteristically related to the transverse electric field via 
cE± = A. This allows us to write Eq. (|2.2I) in a form which treats the transverse and longitudinal electronic field 
on the same footing: 
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E|| •dr+ / E_L-dr) . (2.3) 
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Thus one may view the transverse EM modes as the cause of additional voltage fluctuations V{t) — Ej^ • dr with 
zero mean {V{t)) — 0. This interpretation is physically appealing but still insufhcient to understand in depth the 
detailed nature of the effect of Ej^ on (f>. 

For one thing, we may wonder whether the fluctuating transverse voltage may generate a "slow" contribution to the 
r.h.s. of H2.2|) which would relate to the "fast" part through an appropriate expression of the fluctuation-dissipation 
theorem. This would be analogous to the already noted relation between the second and third terms on the r.h.s. of 
Eq. H2.1|l . In this regard we note that the slow term proportional to cf) in 12.111 is but the Markovian limit of a more 
general retarded expression^, 

rt 

Tqpit - t')<j){t')dt' — > {h/Ae'^R)4){t) . (2.4) 

Clearly, the fluctuating current I{t) plays a role for 4> analogous to that which the fluctuating transverse potential 
V{t) represents for N . This observation suggests that N and (f) satisfy formally similar retarded equations of motion, 
which must be the form 



N{t) ^ -^sinc^{t)~ I Tqp{t-t')4>{t')dt' + -I{t) (2.5) 



m - f rEMit-t')N{t')dt' + ^Vit) . (2.6) 



A rigorous derivation and ultimate justification of this result is possible after the analysis given in Sees. IIIII IIVI and 
fvl where particle-bath couplings are derived which yield Eqs. (j2.5(l and H2.6|l as the Heisenberg equations of motion 
for N and 0. The dissipation kernels rEM(i) and rqp(t) are given in terms of the spectral functions JEM(t^) and 
Jqp(a;), respectively, through the relations 



TEM{t)= JEM(w)cos(wt)dln(tj) (2.7) 
Jo 

/■oo 

rqp(t) = / Jqp(w)cos(tJt)dln(w) . (2.8) 
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The current and voltage fluctuations are related to the spectral functions through the fluctuation-dissipation theorem, 

{vmm - ^ r jemH . '"p^;""!^ (2.9) 

°G J-oo 1 — exp (— ft/5cjj 



{mm) = 2e^ / J^^{cu)^^^^;^pi^du; . (2.10) 



exp^iujt) 
1 — exp {—hfiuj) ' 

As expressed in Eqs. H2.5(l and l|2.6|l . the symmetry between and /(i) on the one hand and N and V{t) on the other 
hand, is most evident. Both stochastic "sources" obey two independent fluctuation-dissipation theorems. The study 
of this type of "mixed" dissipation has been pioneered by LeggettS? for quantum systems. Dissipation due to coupling 
to the momentum variable has been addressed for an ohmic bath in Rcf. The formal symmetry between (f> and N 
is broken not only by the different potential terms, already present in Eqs. H2.5|l and (|2.6(l . but also by the different 
form of the spectral functions Jem(w) and Jqp(ti;) describing the noise. Jqp(aj) is usually taken to be ohmio^Si^ 

^.P(^) = (2.11) 

which through Eq. H2.4|l makes Eq. (|2.5I) become Eq. H2.1(l . It also guarantees that Eq. (|2.10() achieves the correct 
high temperature limit in the form of the Einstein relation mentioned above. By contrast, we expect that, in three- 
dimensional free space, the spectral function of the electromagnetic field is cubic. Actually, a detailed calculation 
performed in Sec. IIIII vields 

J^m{uj) = 3^^, , (2.12) 

where d is the distance between the two electronic clouds and a ~ 1/137 is the fine structure constant. Moreover, the 
spectral function Jem('j-') h^^s to be regularized by an upper cutoff frequency ojc which is the inverse response time 
of the electromagnetic field. This point is discussed in detail in Sec. IIIII where the cutoff frequency is shown to be 
determined by the surface area of the electrodes. Typical values for the diameter of the electrodes of Josephson point 
contacts are about 100 nm. This translates to a frequen cy of lUc — 27r x 10^^ Hz which is higher than typical values 
of the Josephson plasmon frequency tjjp = \/EjEc — \/2eIc/'hC ~ 27r x lO^*' Hz by several orders of magnitude. 
Being the cutoff frequency high, we may approximate Eq. (|2.6|) by a form which is local in time: 

m = -fm + -^m) -^Nit) + -v^t) . (2.13) 

Eliminating (j) in favor of N and viceversa, Eqs. (|2.1|) and H2.13|l become two decoupled fourth-order differential equa- 
tions whose inhomogeneous, fast fluctuating parts are however related to each other. Speciflcally, if we approximate 
sin (i ~ (i, we obtain 



had? ■■■■ , , 4ad^ / Ej hwc \ • • . x / %a(i^Eji^(\ ,-v/ n 1 ,v. x 9 ,w x 

3;^ ^ - ^ (f - 2^ j ^ + ) ^^'^ + ^^^^^ + "'^^^^^ 



U{t) ^V{t) - ^V{t) , (2.14) 



2e ^ ' ^ ' 2eR 

fiad^ 4ad^ [ Ej tuvc \ "j f Sad^EjLOc\ 1 , 2 



3e2Ec2 ^ ^ ' ^c^ \h 2Tre^R J ' ' y Zirfic'^ J ' RC 

En ^, X Sad^Wc ■■, x 2a<P ■■■ , , 2e x ,x 

^^m + -^m - ^ / w + ^vit) . (2.15) 

Here we stress that I(t) and V(i) are two uncorrelated stochastic processes. Nevertheless, their simultaneous coupling 
to the junction generates interference terms. The most obvious signature of such interference is the fourth order time 
derivative in Eqs. H2.14|l and (|2.15|l . since that term vanishes whenever any of the two coupling constants R^^ or a 
vanishes. 

We observe that the left hand sides of Eqs. (|2.14() and (|2.15|l are identical except for an interchange of N and 0. 
This means that the stochastic differential equations satisfled by iV and only differ in the fast fluctuating part. The 
structure of Eqs. (|2.14(l and H2.15|l suggests the introduction of a generalized "response" function 

- \ I \ 2 2 E J f ITT UJc\ Ec ITT / X 

X (t^)=t^jp-t^ +^ ( y-— I ■^EM(w) + -^YJqp(t^) + YJqp(w)JEM(w) ( y-— I . (2.16) 
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A more general form of the response function H2.16|l . valid for arbitrary spectral densities, will be given in Sec. IV Bl 
In Fourier space, Eqs. (|2.14|) and H2.15|l read 



N{uj) = xH 



iLU , ^ 2e / E J in , ^ , , , 



iLu2e 



Ec 
h 



to 
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Ultimately we are interested in autocorrelation functions of the form 

mm')) = ^ / (i</>(c^)p)e'"(*-*'M.. , 



(2.17) 
(2.18) 

(2.19) 



where interference terms must clearly play a role. 

QED effects can be studied more clearly in Eqs. (|2.14|) and (|2.15f) if one neglects the effect of the electronic 
environment, i.e. if one eliminates in Eqs. H2.14|l and (|2.15|) all terms depending on or I{t). One obtains 



AEjacP ■ 



+ ^^)m)+^l.Nit) = -^vit) 

i+^^^)0(o+-,?p^(o - 'ivit), 



(2.20) 
(2.21) 



which are equations of the Abraham-Lorentz type^^. The fact that the EM field couples to the particle number and 
not (directly) to the phase is reflected in the different r.h.s. of Eqs. (|2.20() and (|2.21|l . which are proportional to 
V{t) and V{t), respectively. This means that, although their slow parts obey identical equations, (p and N experience 
different fluctuations terms and thus possess different power spectra. Close inspection reveals that (12 .201) and (|2.21|l 
are the classical equations of motion for a dipole in a harmonic potential interacting with its own radiation. Indeed it 
could have been anticipated that the equation of motion for N should be that of a radiating dipole, since the relative 
particle number A'' couples to the EM vacuum through the electric dipole which it generates, P = 2eNd. In effect, 
neglecting the fluctuating term as well as the mass renormalization terms, Eq. H2.2()|l can be written as 



Pit) 



Pit) - r P (t) = 



(2.22) 



with T — AEjaCp/Zhc^. For a radiating harmonic dipole of charge e and mass m the characteristic time r is given 
b}ii^ T = 2e'^ /imc' . We conclude that the equation of motion for A^ is that of a radiating dipole with effective mass 
TO = h^/3Ej(f. 



III. COUPLING TO THE ELECTROMAGNETIC ENVIRONMENT 



A general description of the electrodynamics in a Josephson junction with electrodes in the x — y plane is given 
by a two dimensional sine-Gordon Hamiltonian for the position dependent phase2& <^(a;, y). The characteristic length 
scale of the theory is the Josephson penetration depth Aj. A small Josephson junction is characterized by an area A 
of the electrodes which is much smaller than Aj. In this limit spatial variations of the phase in the x — y plane are 
energetically highly unfavorable. They are suppressed like ^ ^fA/Xj. Then, the sine-Gordon model reduces to the 
Hamiltonian of a pendulum, first discussed by Anderson^i 

= Ej il - cos4>) + -^N^ + Hem ■ (3.1) 

The phase difference across the junction and the number of transferred Cooper pairs are canonical conjugates: 
\(f>,N\- = i. As indicated in Sec. m we introduce the coupling to the electromagnetic environment replacing (f) 
by its gauge-invariant expression cf) — cjyi — (j)2 + A12, where A12 = (2e/?ic) dr • A(r). The line integral connects 
two points 1 and 2 lying deep enough in the left and right superconductor, where phase fluctuations are completely 
suppressed. That means, it extends over a minimal length of 2Al, where Al is the penetration depth of the magnetic 
field. In dissipative quantum systems one usually assumes the influence of the system onto the bath to be small. 
However in our case the system is made by superconductors which of course alter drastically the low frequency be- 
havior of the electromagnetic vacuum. The screening of the EM field out of the superconductor should be taken into 
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account by an appropriate choice of the wave functions A^a (r) of the normal modes kA in terms of which the vector 
potential A(r) is expanded. These must be the modes which diagonalisc the EM Hamiltonian. Then we can write 



= £A(AA(r)akA + H.c.) , 

kA V J 

Hem = ^ Tit^fcakA^kA ■ (3.2) 



kA 



The influence of the superconductor on the EM field is encoded in the set of orthonormal functions Ak\{r). They are 
solutions of a Helmholtz equation with boundary conditions imposed by the geometry of the junction and the nature 
of the superconducting state. H^m is expressed in its usual form. It is also convenient to expand 

AkA = ^ AkA(r)£A-dr (3.3) 
Shifting (j) by — A12 via a unitary transformation H — > U~^HU with U = exp(iA'^Ai2) one obtains 

H = Ho - NY.i^T^e'nuJkf'^ (zAkAflkA - ^KxA.)) + S^e^TV^ ^ jAkAp • (3.4) 

kA kA 

If the superconductors are equal on both sides one has a reflection symmetry with respect to the x — y plane. This 
symmetry and the trivial observation that A12 = — A21 are sufficient to conclude that all AkA are antisymmetric in 
Z2 — z\ and purely real. This allows us to further simplify Eq. H3.4|) . Finally, with OkA i^kA we arrive at 

H = (1 - cos 0) + ( + ^ ^ J iV2 + TV ^ MkA (okA + a^) + J2 ^^fe^kA^kA , (3.5) 

V kA / kA kA 

where /^kA = (87re^?iijjfe)^/^AkA- The Hamiltonian H3.5() is of the Caldeira-Leggett type with a periodic potential. 
However, there is a subtle but important difference. If, as is standard, we interpret (j) as the position operator and 
N as the momentum, then we must conclude that the EM field couples to the momentum operator instead of to the 
position. This type of coupling is often referred to as anomalous coupling^^^ 



A. QED spectral density 



In order to obtain the spectral density Jem('^)j defined as 



Jem(w) = i^/h^) ^ fJ^x^iuJ - ujk) 



(3.6) 



kA 



one must in principle solve for the normal modes ^kA(r), which are the solutions of a complicated boundary value 
problem. For simplicity, we approximate the EM normal modes as those of free space with appropriate cutoffs which 
account for the real geometry. We expect this to be a good approximation for point-contact Josephson links while it 
will provide only a qualitative description in the case of parallel plates separated by a dielectric. To obtain J^mi^) we 
substitute the sum by an integral, J2kx — * {V/8tt^) J^x J k^dk d{cos9k) d(j)k and take advantage of the transversality 
of the electromagnetic field e\ ■ ek — 0. We choose ei — eg and £2 = e^. With ■ — 0,ee • = —sm9k and 
Gfc • 62 = cos 9k we arrive at 



JEM(t^) = ^ ^ dy [ ^ - l]sm [ — y 



1 



Lod 



1^ 



tod g. / 2Lud\ ( '^^d\ ,■ f '^'-^d\ 

c \ c J V c y Aojd \ c J 



(3.7) 
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>^em(w) has the cubic infrared behavior predicted by the electric dipole approximation^^. In the ultraviolet regime it 
diverges quadratically. Specifically, we have 

>/em(w) = ^ for — >1. (3.8) 

c c 

Despite being a sum of oscillating functions, J(w) is monotonically increasing. The crossover from cubic to quadratic 
behavior takes place on the scale of c/d. Finally, as already mentioned above, the reduction of the two dimensional 
sine-Gordon model to Eq. H3.1|l requires a constant A||. This is taken into account by an ultraviolet cutoff function 
g{u}/ijJc)- Any function g with g{Q) = 1 whose modulus vanishes at infinity faster than any power-law decay serves as 
a cutoff function. One can choose also an algebraic cutoff function provided that limcj-,oo Jem(ijj) = 0. We take the 
cutoff frequency ujc = c/^/A, where A is the surface of the electrodes. For point contact structures one usually has 
vC? > d. Therefore, in order to obtain rough quantitative estimates one may use the simple form of Jem(w) in the 
infrared limit given in Eq. (|3.8|) multiplied by the cutoff function. We express all lengths of the structure relative to 
the cutoff frequency Wc ~ c/a/A. This leads to the introduction of the aspect ratio / = d/^/A. In a point-contact 
experiment, we usually have / ^ 1. Now the spectral density reads 

f UJ 

where the dimensionless coupling constant 

7=-^ (3.10) 

has been introduced, which hereafter will be treated as independent of the cutoff. Eq. (|3.9|l is the spectral function 
we will use in the remainder of this work. Strictly speaking the cutoff ujc — c/^fA only applies to k||. In the direction 
perpendicular to the junction the cutoff length is much smaller and of the order of the Thomas-Fermi screening length 
Atf — lOA. However, since the perpendicular component of the mode wave vector does not couple to 0, we may assign 
to J Akz the cutoff which we wish. In particular, we choose an isotropic cutoff for convenience. 



UJc 



Jem{^^) =1 ^g\ — ] , (3.9) 



B. Equilibrium correlation functions 

The Hamiltonian (|3.5I) is still difhcult to handle due to the nonlinearity of the cos </) term. It has been studied 
extensively in its different aspects and with a variety of different methods^^ . In our present study we are interested in 
the equilibrium dynamics when fcsT, Ec <SC Ej. Then {(jp') <C 1 and we can reasonably approximate 1 — cos (j) ~ ?!'^/2. 
The problem reduces to UUersmas model for a damped harmonic oscillatOB^^i^. For thermal initial conditions the 

symmetrized autocorrelation function of the particle number c'^{t) = ^{[N{t)N{0)]^) and of the phase difference 

C'^\t) = i([0(i)0(O)] + ) can be expressed in terms of the susceptibility^i x(aj) — x' + ix" ■ 

4n^W = ^ y'%"(c.)cos(c.i)coth(^^) 

X~\uj) = - _ ^Ej/n)J^M{^) , (3.11) 
In the last line of Eq. (|3.11|) we have introduced the function 

f{uj) = UJ dw — 

Jo uj' [uj' — w^j ^ 

which is defined for any sufficiently well behaved function f{uj). We note that /(w) is always symmetric, although 
/(w) is not necessarily so. One can consider /(w) as a kind of symmetrized Riemann transform of f{uj). 
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To evaluate the quantities in Eq. H3.11|l we employ the cubic spectral density given in Eq. 1)3. 9|l . which scales with 
u)~'^ . This property guarantees finite results in the Markovian limit ~^ oo. By definition, the cutoff frequency is 
much larger than Ej and the Josephson-plasmon frequency wjp. Then, to leading order in oj"^ the contribution of 
the principal value is just NcLOc, where Nc = /g g{x)dx. Below we take g{x) — 0(1 — x), which implies Nc — 1. At 
zero temperature, the integrals in Eq. (|3.11() can be expressed in terms of the roots of the equation 



n 



2h 



■CJ^tJjp = 



The polynomial has one purely imaginary root iX and a pair of roots iz,iz*. We obtain 



4n^(0) 



1 



TT^-f (A-z) (A-z*) 



Re In 



-In 



< + y 

A2 



Re {z - A) 
Im (z) 



Im In 



UJ^+Z^ 



n^E^jj (A-z)(A-z*) 



Re A-" In 



A2 



In 



UJ^ + z2 



Re(z- A)^ 
Im (z) 
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One calculates the parameters A and z = zi + iz2 perturbatively in powers of Ej/ujc and ivjp/ujc 
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The Josephson plasmon frequency acquires an imaginary part (zi) only to second order in uj~^. On the other hand, 
Z2 — — Re(iz) decreases by an amount EjjujjpNc/2?iuJc, which may be viewed as a macroscopic Lamb shift. Unlike 
in the ohmic case, a crossover to overdamped oscillations is not possible if remains large enough. Expanding 
Eq. H3.15|l in powers of ui~^ one obtains 
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(3.17) 



In the limit 7 — > the susceptibility of the unperturbed harmonic oscillator is recovered and we obtain the well-known 
results for (iV^) and {(f>^). The sign of the correction reflects that the QED field tends to measure N and the expense 
of increasing the uncertainty in 0. 

For finite temperatures it is convenient to introduce an algebraic cutoff function. We choose g{x) = (1 + a;^)~^. 
With this choice the number of zeros of x~^(a;) is finite and we are in a position to evaluate both integrals in Eq. (|3.11|) 
exactly. We obtain the following low temperature dependence. 
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(3.18) 



Note that for lUc 00 there are no finite T corrections. A sketch of the derivation of these results is given in Appendix 

El 

The results in Eqs. H3.17|l and (|3.18|) deserve some discussion. First we concentrate on the mean square of the phase 
difference. One can write the leading contribution to ((/)2) in Eq. (|3.17|) in the following way 



1 Ec_ 
2V Ej 



2a 
3^ 



(3.19) 



This equation is one of the central results of this work. It clearly distinguishes between the uncertainty of the phase 
difference due to the longitudinal electromagnetic field and that due to the transverse electromagnetic field. The QED 



correction depends on the details of the junction only through the aspect ratio /. Junction with / close to unity 
are possible in the case of point contacts. However, the universal constant 2a/?>TT ~ 1.5 x 10^'^ renders the influence 
of the transverse field a tiny effect. The contribution of the longitudinal field contains several quantities which are 
experimentally relevant. First, we have the critical current in Ej — 1ilc/2e. It may be approximated by the formula 
of Ambegaokar and BaratoflSi 
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(3.20) 



at zero temperature. It may be viewed as decreased by the QED field through and effective Debye- Waller factor. 
Second, we have the capacitance in Ec = 4e^/C. 

The case of a plate capacitance is characterized by a/A ^ d > fc^^. This limit corresponds to an aspect ratio / <C 1. 
The capacitance can be approximated by C = eA/Aird. The normal state resistance of a tunneling junction is given 
by = e^N {EF)AkWT\'^ /2, where N{Ef) is the density of states per spin and jTp is the average transmission 



probability for electrons at the Fermi surface, 
approximate 



For a rectangular barrier of height not much greater than Ep, one may 



\T\ 



-bkpd 



(3.21) 



where a and b are numbers of order unity and d' < o? is the barrier thickness. For qualitative estimates we may 
assume d! — d which underestimates slightly the QED effects. For simplicity, we also assume a ~ 5 « 1 in the ensuing 
discussion. If we only include the leading QED correction, Eq. H3.17|l becomes 
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(3.22) 



where dp = kpd, Ap = Ak'^p, uq is the Bohr radius, and = hvp/nA is the zero temperature coherence length. Since 

typically ^ 2000 oq, the ratio between the two terms within square brackets is {150 / a) d ^ ^ ex.p {dp / 2), which 
takes a minimum value of ~ 10^ at dp — 3. We conclude that, for plate capacitances and within the approximation 
of free space photons, EM vacuum fluctuations contribute negligibly to the quantum phase spread. 

The situation is not much better if one considers the case d \/A, as might correspond to e.g. a tunneling point 
contact. The capacitance is then approximated as C ~ \/A. Eq. H3.22|l is replaced by 



3nAi 
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(3.23) 



For the QED term to be comparable to the charging contribution, one would need unrealistically small areas which 
would effectively suppress the supercurrent. 

The case of a planar junction which has been considered in the first place presents the unfortunate property that 
both the longitudinal and transverse contributions scale identically with the junction area. That makes it difficult 
to overturn the natural smallness of the QED correction. We note however that we have considered the spectrum 
of photons in free space, an approximation that was admitted to be less adequate for the case of broad Josephson 
junctions. In fact, the structure of Eq. H3.11|l suggests the possibility of designing devices where the coupling of the 
phase to the photon field could be amplified. 



IV. COUPLING TO THE ELECTRONIC BATH 



A. Effective Hamiltonian 



A complete treatment of the macroscopic phase dynamics in Josephson junctions must include the effect of quasi- 
particles and Cooper pairs. It was shown in Ref. ^ within a path-integral language that the electronic degrees of 
freedom can be taken into account by a coupling to two independent oscillator baths. Starting from a microscopic 
tunneling Hamiltonian, they derived an effective Lagragian in which the phase difference couples to two heat baths 
which we label with superindices (±). Here we propose an alternative Hamiltonian-equation-of-motion approach to 
derive the coupling and spectra of these effective baths: 
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(4.1) 



The coupling constants A^^^* are defined tlirough the spectral densities they must yield. Specifically, the expression 



(4.2) 



must be equal to 



where 



tanh 



Nqp{uj) = N{Ef)Q [huj - A) 
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(4.3) 



(4.4) 



is the quasiparticle density of states^SiSS. 

We derive this Hamiltonian in subsection IIV B| avoiding the path-integral formalism used in Ref. Our criterion 
will be that the heat baths must yield the correct equations of motion for the phase 4>{t). The reader who is familiar 
with the subject or is not interested in this point might skip the next subsection an take the Hamiltonian H4.1|l for 
given. 



B. Derivation of the effective Hamiltonian including coupling to electronic degrees of freedom 



Our starting point is the standard tunneling Hamiltonian for two coupled superconductors in the Bogoliubov-de 
Gennes mean-field approximation22i22i2L2SiS. 



BdG — 
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(4.5) 



with the 4x4 matrix 
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-T*^ r' 
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(4.6) 



The order parameters A(r, r') obeys the self-consistency condition within the bulk of each superconductor. In Eq. H4.5() 
we suppressed the dynamics of the electromagnetic field for momentary convenience. The operators H{r, r') as well 
as A(r, r') — |A(r, r')| exp(i0(r, r')) are local in space. After a further gauge transformation 



C/ = exp|^ I dr*t(r)diag [</)(ri), 0(ri), (/.(r2), (/.(r2)] *(r) 



and in the grand canonical ensemble, the entries of ft in Eq. H4.6I) look explicitly as follows: 



Ai^fl(r,r') = |A«,i.(r)|(5(r-r') 

Tr.^r^ = exp(z(/.(ri,r2)/2)fr,,r2 



5{v - r') 



(4.7) 

(4.8) 
(4.9) 
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We remember here that A'' is the semi-difFerence of Cooper-pairs between the left and the right electrodes, N = 
(TVi — N2)/2. We have written the tunneling part of the Hamiltonian Ht in such a way that the gauge dependence 
becomes explicit. 

Tr^r' f;,r'expiA(r,r') 

A(r,r') = ^ /" A(r)-dr= -A(r',r). (4.10) 
he J^, 

Try is gauge invariant and, without loss of generality, can be assumed to be real symmetric. Thus the fermion-fields 
Tpa are the only quantities in Eq. (|4.5|) which break local gauge invariance. The notation above may look unnecessary, 
since in the standard descriptions of Joscphson junctions'^'' the tunneling is treated locally and thus the gauge term 
in (|4.1U|) can be safely dropped. However for our purposes it is convenient to keep the gauge invariant line integral 
in the definition of the tunnel matrix element. The virtue of the transformation H4.7|l is that it permits to elucidate 
the role of the phase of the order parameter, which appears in the transformed Hamiltonian only in a gauge invariant 
form, through either v(r) or 0(r, r') = ^(r) — (/'(r') + A(r, r'). 

In order to derive an effective theory for the phase one makes two further approximations. First, one assumes 
for simplicity that the tunneling matrix element is independent of the electron momenta, i.e. one sets r(r,r') = 
r^(r|| — r'|)5(2;)^(z')iMi. Second, one approximates the unperturbed parts i?L^fl(r) substituting v(r) and A(r) by 
their mean field bulk values v(r) = 0, z ^ 0, A(r) = A. This is a reasonable approximation in the tunneling 
]^jjjjj^23j37j38 Using the above approximations we may write the tunneling Hamiltonian as 

H = {Ec/2)N^ + Hl+Hr + Ht 

i(r),(T ;(r) 

Ht = rcos(0/2) (cLcra + H.C.) + zrsin((/)/2) ^ (c^c,, - H.c.) . (4.11) 
In the third line we have used the definition of the Bogoliubov quasiparticle operators ^k,^ = UkCk\ — VkC^_f.^ and 

ik,i = ukc-ki + Vkcikr 

Now we adopt a bosonization procedure whereby the Hamiltonian H4.11|l is mapped onto an effective Hamiltonian 
of the form 

H.. - + cos(0/2) Y + + ,i,(^/2) J- ^t' it' + t'') + E ^^^''t' ' (4-12) 



where the fe^^-* are bosonic operators and the coupling strengths a'^'' are yet undefined. To further specify the 
Hamiltonian H4.12|l we use the following strategy: We derive the Heisenberg equations for the canonical conjugate 
variables N and as well as for the electronic degrees of freedom of the Hamiltonian H4.11|l . Then the equations 
for the electronic degrees of freedom are solved in perturbation theory. Upon elimination of the electronic degrees 
of freedom one obtains two coupled integro-differential quantum Langevin equations'^^ for the Heisenberg operators 
N and (j). It turns out that, up to second order in the hopping T, these equations have the same structure as those 
obtained from Eq. (|4.12|) . This allows us to adjust the coupling parameters a|^^ in Eq. (|4.12|) in such a way that the 
equations of motion for N and (f> become identical. Below we implement this strategy. 
From the Hamiltonian (|4.11|l we obtain the Heisenberg equations 

<j){t) = h-^EcN{t) 

^(^) = In ^ cL(i)c..(Oe^^(*)/^ + H.C. . (4.13) 

To close these equations we evaluate the time evolution of the electronic part perturbatively in Ht- To first order we 
obtain 



J — oo 1 



12 
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(4.14) 



and proceed similarly with the other relevant fermionic operators. The superscript / denotes operators in the inter- 
action picture. Introducing Eq. (|4.14|) into Eq. I|4.13f) yields 



dt'a{t - t') exp i-^-^ 



' dt'/3(t-t')exp('zM±^^)+H.c. 



(4.15) 



Since the uncoupled Hamiltonians Hl j^ in Eq. H4.11|l are quadratic forms, the corresponding grand canonical ensemble 
is a Gaussian ensemble. Therefore the effect of the operator 



Ht) = ^Y.^\l{t)cUt) 
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(4.16) 



l.r.a 



is determined by its second moment {F{t)F{t'))B, where subscript B denotes average over bath coordinates. These 
quantities are related to the thermal averages of the operator valued dissipation kernels 

IT 
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(4.17) 



by a fluctuation-dissipation theorem. In order to derive a differential equation for the the iV-point function 
{(j>{ti) . . . (f>{t]sf))B from the Heisenberg equations Eq. H4.15(l . one must in principle evaluate averages of the form 
{a(ti) . . . d(ijv))B- However, since all averages are made over a Gaussian ensemble. Wick's theorem applies and they 
reduce to products of {a(ti))B and (/3(ii))B respectively. Therefore, in order to completely determine the quantum 
Langevin equation Eq. (|4.15|) . it is sufficient to calculate the averages (d(t))B and {(3{t))B- These quantities are 
calculated in standard textbooks on superconductivity^^-'^*'. The outcome is the well known result that {a{t))B and 
0{t))B are the sine-transforms of the sum and difference, respectively, of the two spectral functions J^^^{uj), 



{a{t))B = a{t) = 

CmB^m = 



J^+\uj) - /^Hw)] sin(wi)da; 



(4.18) 



where J'^^\ijj) are the spectral functions defined in Eq. (|4.3|l . 

We can also derive a quantum Langevin equation for N and (j) from the bosonized Hamiltonian Eq. (|4.12(l . Since 
the interaction here is linear in the bosonic operators, the procedure is exact (by contrast, fermionic correlators 
are computed up to second order in Ht)- The resulting Heisenberg equation is identical to that in Eq. (|4.15|> ^S^. 
However, in this case the response kernels are given by 



i 



A 



i+V 



A 



(-)' 



(4.19) 



Since the coupling parameters aI^"* have not yet been specified we are free to choose them in a way that adjusts the 
dissipation kernels in Eqs. (|4.19|) to the kernels (|4.18l) derived from the microscopic theory. Thus we may conclude 
that the two Hamiltonians H4.11|l and 1)4. 12(1 . together with the spectral density H4.3|l are equivalent to order 0{T^), 
at least as far as the dynamics of the phase is concerned. We stress that, due to the Gaussian nature of the averages, 
this equivalence holds not only for the dynamics of the mean (0(t)) and {N{t)) but also for the dynamics of the higher 
correlators. 
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In order to make contact with the Hamiltonian Eq. (|4.1() we add and subtract in Eq. H4.12|l the expression 



i ^ i ^ 

2 ^ cos + - ^ , (4.20) 

i i 

which allows us to complete the square. The </> dependent term in Eq. H4.20|l can be calculated using the definitions 
of A^^'' in terms of the spectral densities J*-^^ given in Eq. (|4.3|) . We obtain the result of Ambegaokar and BaratoflSi, 
namely, 

1 A(+)' - A(-)' , 1 r ^ - j(-) , 

- > cosm = - duj cosffl 



C0S( 



= ^tanhf^) 
4e2i? \ 2 J 

= E J cos (p. (4.21) 

The second term in the r.h.s. of Eq. H4.20|l is divergent but, fortunately, it is independent of 4>. Since we are allowed 
to add an arbitrary constant to a Hamiltonian, we subtract this infinity and keep, for convenience, only a constant 
Ej from this second term. This yields the Hamiltonian (|4.1|l . 

Of course the two Hamiltonians (|4.1|) and H4.12|l are identical, but the physics is more transparent in Eq. H4.1|) . It 
is the Hamiltonian of a particle moving in a cosine potential coupled nonlinearly to two independent bosonic baths. 
Since the interaction is in the form of a complete square, the coupling to the heat bath does not change the potential 
landscape. This is a well known fact in the theory of dissipative quantum systems^i^. 



V. SIMULTANEOUS COUPLING TO PHOTONS AND QUASIPARTICLES 
A. Mixed coupling Hamiltonian 



If we now introduce again the dynamics of the electromagnetic field as explained in Sec. IIIII we arrive at a rather 
comprehensive description of the dynamics of the phase in a Josephson-junction. It is instructive to cast by a further 
unitary transformation the resulting Hamiltonian into a highly suggestive form. Applying HU with 



U = exp 



cos 



1 1: 
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(5.1) 



onto H in Eq. (|4.1|l shifts the interaction of the "electronic" bosons with (f) into an apparently more complicated 
interaction with N. However, we can write the transformed Hamiltonian as 



H = ^{N + SNf + Ej[l-cos{(f) + S(j))] 



(5.2) 



kA 



This form is rather general and physically quite elucidating. We observe that the position degree of freedom of a one 
body Hamiltonian acquires a fluctuating part. So docs its canonical momentum. These fluctuating parts obtain their 
own dynamics as bosonic excitations. The specific form of these fluctuating operators has either to be derived from a 
many-body Hamiltonian or it has to be modelled phenomcnologically. In our case the fluctuating term of the particle 
number 



SN ^ - sin 
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(5.3) 



is due to the hopping of Cooper pairs and quasiparticles across the junction. It is phase dependent and was derived 
from a microscopic model. An interesting message from Eq. (|5.2|l is that the Coulomb interaction does not distinguish 
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between N and SN. It is only sensitive to the overall particle number difference. This may be viewed as a form of 
minimal coupling, since it reflects the "gauge invariance" of the Coulomb term and is analogous to the criterion 
whereby the coupling to the EM field is introduced by replacing (f) by its gauge invariant form. The fluctuating part 
of the phase is Scj) = A12 [see Eq. 

Despite the appealing compactness of Eq. (|5.2|l . for calculational purposes it is more convenient to shift both N 
and (j) by their fluctuating part ~6N and —dcf) through a unitary transformation. We employ the inverse of (|5.1(l to 
shift N and U ~ exp(— «iVAi2) to shift (j>. The resulting Hamiltonian is much closer to the standard Caldeira-Leggett 
form 
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«kA + ^iV),(5.4) 



with the coupling constants A^ defined by J'^^'(w) through Eqs. (|4.2|) and (|4.3|) and /ikA defined by J-em{^) 
through Eqs. (|3.6|l and (|3.9|l . The nonlinear coupling of the phase to two different baths is rather difficult to handle. 
Following Ref.^'^, one considerable simplification is achieved by keeping only a local contribution in the (3 kernel 
integral of Eq. (|4.15|) . This is a valid approximation if the phase varies slowly on the time scale in which decays (3{t), 
which is set by 7i/A. Yet another simplification is obtained if one keeps only the first non-vanishing term of the 
exponential in the a kernel integral in Eq. I|4.15|l . With these approximations the nonlinear coupling of the phase to 
the electronic degrees of freedom has been replaced by a linear coupling to only one effective heat -bath of an "almost" 
ohmic character. These are the basic assumptions which underlie the RCSJ model. The limits of validity of these 
approximations have been discussed in the literatureSi^. In the RCSJ limit, the Hamiltonian H5.4|l becomes 

Ec 
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The coupling parameters Ai are adjusted to satisfy Jqp(ct;) — 2h 
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— tanh 



(5.5) 



T (-■-!)]} 



(5.6) 



K being the normal resistance of the junction. This means that, in the RCSJ approximation, diffusive effects due to 
/3(t) are suppressed. 

In the high temperature limit, Eq. (|5.6() yields the ohmic spectral density for a normal conductor junction given 
in Eq. (|2.11|1 : Jqp(w) ~ ^/Rq, Rq = ^TiRe^ /%. This ultimately justifies the RSCJ model. A plot of Jqp((^) for 
various temperatures has been given by Harris^. One observes that, already for relatively high temperatures, the 
spectral density deviates from strictly Ohmic behavior. In particular for ksT ^ A/Ti the spectral density is better 
approximated by the zero temperature spectral density, which can be calculated analytically*^, the most prominent 
feature being the existence of a gap at w = 2 A. It is known that an infrared gap in the spectral function reduces the 
dissipative effect of the bath. Therefore, in the RCSJ model the effect of quasiparticle tunneling is overestimated. 

A central physical feature is that the two baths in Eq. H5.5|) have competing effects. On the one hand, the 
quasiparticle environment tends to "measure" 4> (the position of the particle). On the other hand the environment 
"measures" A^ (the momentum of the particle). Of course, due to the uncertainty principle, the two baths cannot 
be equally effective. In general, it is not clear a priori which mechanism will dominate. In this case, however, both 
because QED is weak and supcrohmic, one expects the quasiparticle measurement of the phase to be dominant. 

Apart from its relevance in the present context, the model which we have constructed is worthwhile to be studied 
in its own right^^. We address this subject in a more general context in Ref4&. Here, however, we restrict ourselves 
to the case in which one bath is ohmic (quasiparticles) and the other one is cubic (photons). 



B. Equilibrium correlations 



Like in Sec. IIII Bl we use the harmonic oscillator approximation for the nonlinear potential terms. In contrast to 
the case of a harmonic oscillator coupled to a single bath, this model is not exactly solvable for an arbitrary initial 
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condition. Fortunately, we are mostly interested in the equilibrium quantities (t) and CjZ^ (t) 



|x([j)|2 cos{Lut) coth{hf3uj/2) 



ImJqp(cj) ( — ^ - 2— ReJEM(w) + \ Je,-m{i^)\'^ ) + c^'^ImJEivilw) \ dw 



|x(cj)P cos{ujt) coth(n/3w/2) 



IuiJemH yji - 2— ReJqp(w) + |Jqp(cj)|^j +CJ^ImJqp(w)| dcj 

Here, we have introduced the generalized "susceptibility" 

X"^(w) = Wjp - - ?i"^£',/JEM('^) - h.~^EcJq,p{'~^) + Jqp(a;) JemC^^) 



(5.7) 



(5.8) 



Again, /(w) denotes the symmetrized Riemann transform as defined in Eq. (13.12(1 . To leading order in to^^, Eq. 
(|5.8|) reproduces Eq. (|2.1t)|) . which was derived from a simplified discussion. Of course, these expressions reduce 



to Eqs. H3.11|l for Jqp(w) = 0. In Eqs. (|5.7|) the spectral function Jem(w) is defined by Eq. H3.7|) . while for the 
quasiparticles we use the ohmic form Jqp(ti;) — {uj / RQ)g{Ljj /VLc) keeping in mind that it overestimates dissipation at 
low temperatures. Moreover we introduce formally an electronic cutoff fic, which we may choose as the conduction 
band width. For the smallest geometries available one has approximately fie ^ ^c- In the following we restrict 
ourselves to T = 0. As pointed out in Sec. lIII Bl this should be a good approximation at sufficiently low temperatures 
ksT < A ~ keTc < n^. 

The evaluation of the integrals in Eq. (|5.7() involves the solution of a fifth order polynomial. Three of its solution 
scale as ~ Wc while the other two solutions shift the Josephson plasmon frequency into the complex plane. Eventually 
the real part disappears completely, leading to overdamped oscillations which are characteristic of ohmic and subohmic 
spectral densities. For {(j)^) we find, to leading order in lo'^ and 
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(5.9) 



where we have used the notation k = {ir / 4Rq) ^/Ec/Ej . Moreover, we have introduced Vandermonde's determinant 
Ajv(r) = Yli^j [xi — Xj), i, j — 1, N. The dimensionless quantities Aq, Ai, A2 are the coefficients of the leading order 



term in lUc of the three solutions of the fifth order polynomial which scale with ujc 
the third order polynomial 



TT 2 y TT TTJ 







They are given by the solution of 



(5.10) 



Remarkably, the roots Ai only depends on the ratio Rq/j- The expression for {(fP) in the first line of Eq. H5.9|l consists 
of two contributions. The first term is recognized as the mean square of the position operator of a harmonic oscillator 
with ohmic damping^!. It is the leading order contribution in fl^ of quasiparticle tunneling. The second term, J^, 
stems from the electromagnetic environment. Obviously the exact analytic form of stated in Eq. H5.9|l is not too 
illuminating. To get further insight we consider the limiting cases Rq/j <^ 1 and Rq/j 3> 1, which allow for an 
expansion of the logarithm in g^j, and a perturbative solution of Eq. (|5.1U|) . For Rq /j ^ 1 we obtain for 



7 
4 




(5.11) 



As we have seen above, the parameter of the coupling to the electromagnetic field 7 = Sa/^/STr is always small. 
However, the dimensionless resistance can be either large or of order unity depending on the experiment. Assuming 
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Rq ^ 1 or, equivalently, k ^ 1, we can also expand the first relation in Eq. H5.9() . We obtain, to first order in Rq^, 



1 [E^ Ec 7 
2\ Ej AEjRq 4 



4Rr 
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TT 

T 



(5.12) 



This result illustrates nicely the competing character of the effects of the quasiparticle bath on the one hand and 
the photonic bath on the other hand. The second term in Eq. H5.12|l is the dominant quasiparticle contribution for 
ilc — > oo. ft leads to a reduction of {(jp) or, equivalently, to an enhancement of coherence. The third term in Eq. (|5.12|) 
is the QED correction to first order in uj~^, which was already found in Eqs. H3.17|l and H;H9|I of Sec. lIII Bi This term 
is always positive. One can estimate the quasiparticle term by using again the Ambegaokar-Baratoff formula (|3.2(J|) 
and the capacitance of a plate capacitor. One finds that the QED contribution would dominate for an unrealistic 
thickness d > 12fic/eA « 100//m. 

For the sake of completeness we also discuss the opposite limit Rq/^ ^ 1. After some algebra, we obtain that the 
leading order term of the "mixing" function is 



7 

Rq 



Expanding Eq. (|5.9|) for Rq <C 1 we obtain 



Rq 



4 In 



2 

R^ 
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The result for {N'^) is formally similar to that of 



1 lEj ln[(K + 1)/(k-Vk2- 1)] 

2 V Ec tta/k^ — 1 
1 . ( ^ 
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(5.13) 



(5.14) 



(5.15) 



2Xf 8Rq 2Rq 
For small EM coupling {j/Rq <C 1) the mixing function /at reduces to 

= 4RI 

We reproduce in Eq. (|5.15|1 the logarithmic dependence on the quasiparticle cutoff frequency, a well known feature 
of the coupling to an Ohmic heat bath'^-. 



(5.16) 



VI. SUMMARY 



We have explored the role of the quantum electrodynamic field as a possible agent of zero-point decoherence in 
some physical systems. We have focussed on the fundamental problem of the dynamics of the phase-number variable 
in a Josephson link. The transverse electromagnetic field has been shown to couple to the relative Cooper number 
through the electric dipole which it generates. The effect of quasiparticles has been approximated with an effective 
oscillator bath which couples to the phase. This poses a new quantum dissipation problem of a type which has 
been little studied: that of a quantum particle coupling through its position and momentum to two different baths. 
The resulting dissipative dynamics of the macroscopic phase-number variable has been investigated by means of a 
Hamiltonian equation-of-motion approach. We have found that the r.m.s. values of the phase and number contain 
interference contributions from the photon and quasiparticle bath. The effect of the QED field has been found to be 
quite small compared with charging effects due to Cooper pair and quasiparticle fluctuations. The unfortunate fact 
that the contribution of the longitudinal and transverse fields to the phase variance scale identically with the junction 
area precludes a simple observation of the effect of electromagnetic vacuum fluctuations in generic Josephson devices. 
However, the sensitivity to the photon and quasiparticle spectra which we have found suggests that the design of 
special structures, where the effect of electrodynamic zero-point fluctuations is amplified, cannot be ruled out. 
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APPENDIX A: DERIVATION OF EQS. ifsl^ 

With the Drude type cutoff g{x) — (1 + the inverse susceptibility involves a fourth order polynomial in the 

numerator and a second order polynomial in the denominator 



(TtEj 2 212 /7\ 2 22 



(Al) 



The roots of the fourth order polynomial come in two pairs of complex numbers iX ~ iXi ± A2 and iz = izi ± Z2. The 
roots of the second order polynomial are given by Aq = Wc exp(±i7r/4). One can evaluate Eqs. I|3.11|l by a contour 
integral similarly to the Drude model^^i^. However, here we will take advantage of the general expressions for the 
mean squares {N'^) and {(jx^) as derivatives of the partition function of the harmonic oscillator^. We have 



{N^) = lnZ(/3) 



- -2i;^(--d^ + ^^^)i^^(^) ' (^2) 

with the partition function of the damped harmonic oscillator 

Z{P)^-^f{^lx{M , (A3) 

We have used the notation Vn = 2Trn/h(3 for the Matsubara frequencies. 

At this point it proves useful to use an algebraic cutoff. That yields a finite number of roots which satisfy a series 
of algebraic relations, called Vieta relationgii. This allows us to express the infinite product in Eq. IjASp compactly 
in terms of F-functions 

?./3o.jP |r(AHp|r(z/z.)p 

4^2 |F(Ao/i^)P ^ ' 

The real and the imaginary parts of the roots A, z scale quite differently with the cutoff frequency. They are approx- 
imately 

Al = l=cj, + 0{lj-^) 

A2 = 7=^ ^ (Ji'^c ) 

V2 M 16^/2^2 ^ " ' 

Using the asymptotic expansion of the F-function we get an expansion of Z{f3) in powers of ksT. 

00 

Z{f3) = + J2 Z^"\P) {ksTT ■ (A6) 

n=l 
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In zeroth order we recover Eqs. H3.17|l with Nc = n/2^/2. The next two orders in the expansion vanish identically due 
to the Vieta relations. The first non-vanishing correction term reads 



1 



360w6u;6p 



+ A*') |z|6 + (z3 + |A|6 - ujIp (Ag + Xf) 



(A7) 



Using Eq. (|A2() and retaining only the highest order term in lOc we finally arrive at the result for {N'^) stated in 
Sec. nil Bl with a quartic finite temperature correction. Moreover we observe that for (0^) the quartic contribution 
vanishes as lOc — *■ oo. Therefore, in order to obtain the first finite temperature correction for (if)'^), we have to expand 
Z{P) up to fifth order in /3. The result with a sixtic finite temperature correction is given in Eq. (|3.18() . 
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